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Abstract. This paper is divided in two parts. In the first part, the inverse spectral 
problem for tight-binding hamiltonians is studied. This problem is shown to have an 
infinite number of solutions for properly chosen energies. The space of such solutions 
is characterized by a hypersurface in the space of hopping amplitudes (i.e. couplings), 
whose dimension is half the number of sites in the array. Low dimensional examples 
for short chains are carefully studied and a table of exactly solvable inverse problems 
is provided in terms of Lie algebraic structures. With the aim of providing a method 
to generate lattice configurations, a set of equations for coupling constants in terms 
of energies is obtained; this is done by means of a new formula for the calculation 
of characteristic polynomials. Two examples with randomly generated spectra are 
studied numerically, leading to peaked distributions of couplings. In the second part 
of the paper, our results are applied to the design of bent waveguides, reproducing 
specific spectra below propagation threshold. As a demonstration, the Dirac and the 
finite oscillator are realized in this way. A few partially isospectral configurations are 
also presented. 
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1. Introduction 

The continuous interest in inverse spectral theory [1, 2] has reached many important 
results connected with cavities, billiards and their inherent quantum chaology [3,4]. 
Apart from the purest interest of Kac’s question and the shape of drums, fresh 
applications of quantum mechanics seem to draw the attention to similar problems 
in different areas such as propagation in crystals. Recent examples in experiments 
connected with the quantum behaviour of matter [5,6] illustrate the important role of 
tight-binding lattices and their shapes. This is also true for other emulations of solids 
in mesoscopic physics, which sacrifice quantumness in favour of geometrical flexibility, 
with the aim of exploring more and more effects inexpensively [7-10]. Along with 
these examples comes the necessity of dealing with inverse spectral problems defined on 
lattices of any dimensionality, in particular tight-binding chains. Thus, inverse lattice 
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design can be conceived as an answer to the problem of a discrete shape producing a 
specific spectral feature. 

The incarnations of tight-binding models are of a great variety, but the focus of this 
paper shall be centered in waveguides, since they can be regarded as elongated variants 
of the more traditional drum. See figure 1. Waveguides can also be open at their ends, 
which makes them accessible through scattering, with slight perturbative effects in their 
resonating frequencies when compared with their closed counterparts. The question is, 
can we hear the shape of a pipe? 

In this paper we provide a general answer to the following two practical questions: Is 
it possible to produce an arbitrary finite spectrum by the mere use of a one-dimensional 
chain without local potentials? Are these results applicable to the transport properties 
of a scattering device? The answers shall be delivered by characterizing the space of 
all tight-binding hamiltonians that produce a suitably defined spectrum, together with 
their plausible realization in bent waveguides operating with a well-established binding 
mechanism at corners [11-14]. 

Similar problems have been proposed in previous works; notably, the system of 
interest has a close resemblance with quantum graphs, for which there are extensive 
studies [15-17]. However, there are also sharp differences between our waveguides and 
one-dimensional wires conforming a graph. For instance, evanescent transport is possible 
only in the presence of a non-vanishing thickness, and it is through this mechanism that 
we shall establish finite difference equations. In this respect, it is important to mention 
that significant results in the field of graphs can be found in [18], relating the spectrum of 
a continuous one-dimensional Laplacian with that of a finite difference operator defined 
on graph vertices. Nevertheless, the design of a specific graph producing an arbitrary 
spectrum is rarely given in the literature. 

We present our discussion as follows. In section 2 we characterize the spaces of 
hamiltonians for a given spectrum and we generate a set of algebraic equations that 
relate hopping amplitudes (couplings) in terms of energies as parameters. This part is 
of a mathematical nature and it is carefully presented. As an application, in section 
3 we study bent waveguide realizations of tight-binding arrays. We demonstrate that 
exactly solvable systems of the inverse problem can be emulated by explicit design of 
the lattice: the Dirac oscillator, the finite oscillator and partially isospectral chains are 
presented. We conclude in section 4. 

2. The problem of inverse design in nearest neighbour models 

We shall work with nearest neighbour tight-binding systems of identical sites, which 
belong to the wider class of Jacobi matrices. We must mention that previous work 
such as [19] and references therein, deals with the inverse spectral theory of such 
operators. Although the theory has more general hypotheses related to all posible 
tridiagonal matrices, the claimed solution in such references rests on the assumption of 
two known spectra for two Hamiltonians belonging to the same sequence - see Corollary 
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Figure 1 . Systems described in this paper: bent waveguides with multiple corners 
and tight-binding arrays of resonators. 


2.5 in [19]. The intertwining of levels and the specific form of the spectral function 
allow to determine univocally the elements of one of the Hamiltonians (either the first 
or second in the sequence). However, our inverse problem is restricted to the knowledge 
of only one spectrum. We shall see that this explains the existence of infinite solutions. 
The determination of such a space of solutions can be seen as a complementary result 
to the known literature. 

Let us establish our mathematical notation: We choose units such that h = 1, we 
denote by N + 1 the number of sites in the array and we assume that any state \<p) can 
be written as a superposition of a complete set of localized (atomic) states \n): 

N+l 

\<P) = ^2<Pn\n). ( 1 ) 

71=1 

We consider the wavefunction (x\n) to be concentrated around the site located at x n 
(e.g. a Wannier function [20,21], but without band index for simplicity). Since we 
deal with arrays of identical sites, we must have identical localized functions related by 
translations ( x\n — 1) = (x — x n + x n _i\n). The atomic states themselves | n) can be 
translated by means of a unitary operator T 


T\n) — \n — 1), T^\n) — \n + 1). 


( 2 ) 
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The matrix elements of T are clearly given by (n\T\m) = d n m _\. We can define an 
observable N associated to the site label n with the property 

N\n) = n\n), (n\N\m) = n8 n ^ m . (3) 

For any function F of the operator N alone, we must have 

F(N)\n) = F(n)\n ) = F n \n). (4) 

The stationary Schrodinger equation corresponding to nearest neighbour chains with 
variable couplings but without local potentials is thus 

H\<f> k ) = E k |0 fc ) = [j F{N)T + T*F(N)] \<j> k ), l<k<N + l, (5) 

where now F(N) encodes the information on the variable couplings, according to their 
location in the chain. In fact, (5) can be written as a recurrence if we set 

N +1 

\<i>k) = '52<i>n\n), (6) 

71 = 1 

from which the following relation holds: 

+ l < k < N + 1, l < n < N + 1. (7) 


Since the system is finite, the previous relation is supplemented with fictitious boundary 
conditions (fro = 0v+2 = 0. Similarly, the couplings to fictitious sites must vanish 
Fo = F n+ i = 0. 

With these conditions, the physical block of the system is determined by the 
non-vanishing elements of the operator F(N). Since we are interested in completely 
connected chains, we impose F n ^ 0 for 1 < n < N, but it is important to stress that 
N can be arbitrarily large. 

A few results can be established already: The spectrum a(H) - with H as in (5)- 
is not quite arbitrary. 

Proposition 1 If E k G cr(H), then —E k € <7 (if), i.e. cr(H) is symmetrical around the 
origin. 

To prove this assertion, we multiply by (—l) n+1 both sides of (7) and define ?/ ; n = 
(—1 ) n 4>ni leading to 

-E^^Fn^ + F^ti, 


( 8 ) 
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which implies that I'lpk) = V ; nl n ) is an eigenstate of H with eigenvalue —E k □. 

A simple corollary of this result is that for N even, 0 € a(H ), i.e. at least one E k 
should vanish. 

Another result of interest here is the nature of wavefunctions in connection with 
couplings. 

Proposition 2 Let H = \F(N)\T + Tt|F(iV)| with |F(7V)||n) = |F n ||n). Then 
a{H) = o{H). 

The proof is given by a gauge transformation when F n is allowed to be complex. If we 
set F n = e lSn \F n \, A n = J2j=i an d multiply both sides of (7) by e* An , we find 

Ek(e iA "<P k n ) = F„(e iA "VS + i) + (9) 

From this relation it follows that the new state 

N+l 

\^k) = elAn( t ) t\ n ) ( 10 ) 

n= 1 

has energy E k with hamiltonian H □. 

Therefore, without loss of generality, we can work with real and positive couplings 
F n in one-dimensional systems without loops. Although we are not concerned now with 
a non-trivial topology of the chain, it is worth mentioning that the influence of magnetic 
fields would destroy (9). 

In passing, we note that in certain applications related to transmission of 
information along chains (e.g. bosonic or fermionic matter in optical lattices [5,6]) 
one needs a second-quantization scheme. The formalism previously developed can be 
generalized easily to such situations: The local-held replacements (f) n (->• d n , 0* (->• a ) n , 
with canonical relations [a n ,atj] = d n m or {a n , a^} = d n m lead to spectral problems 
similar to (7). 

Returning to our problem as stated in (5), we recognize that H can be identified 
with the position operator X of an /- or q- deformed oscillator (e.g. deformations of 
Heisenberg algebras a la Man’ko [22-24]). This type of quantum optical application 
requires an operator N — 1 representing the occupation number of photons at a given 
frequency, i.e. a normal mode in a cavity. What we note from this connection is that 
indeed, position operators have symmetric spectra about the origin and that eigenvalues 
can be distributed either continuously, discretely, in bounded regions, in closed intervals, 
or in the whole real line. 

Finally, the inverse spectral problems that we tackle in this paper must be fed fwith 
the following information: 

i) A symmetric, but otherwise arbitrary set a(H), 


f An exception may arise when one of the eigenfunctions of H is known from the outset. We touch 
upon this case in other sections. 
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Formula (19) 
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Figure 2. Diagram showing the maps between isospectral tridiagonal hamiltonians 
(lower row) and their characterization in terms of X (upper row). 


ii) A set of real coupling constants F n , 

iii) An arbitrary but fixed size N + 1. 

With this in mind, we characterize the spaces of solutions in what follows. 

2.1. The many solutions of the inverse spectral problem for variable couplings 

In this section we determine the space S of all possible nearest-neighbour hamiltonians 
that lead to the same spectrum - complying, of course, with the requirements of 
the previous section. We shall see that such spaces can be characterized by orbits 
generated by linear operators, which however do not correspond to the whole space of 
(. N +1) x (N +1) unitary matrices due to the nearest-neighbour restriction. We shall see 
as well that the gauge transformations defined by (10) are but a trivial part of S. For 
a positive integer n. let us denote by IA (■ n ) the group of unitary n x n matrices and let 
us define Uh(ji) = {UHW\U G U (n)}. Then, for a fixed cr(H) we have an equivalence 
class S h containing H such that 


[Ujf(l)P +1 C S h C Uh(N + 1) 


(U) 


where the contentions are proper except for = 0 or vanishing couplings. The process 
by which we shall work consists in the introduction of two auxiliary operators A and 
X ; A encodes the information of the couplings and A" connects all the hamiltonians 
isospectral to the homogeneous chain. This is shown diagramatically in figure 2 and we 
proceed to explain it. 

First, we determine the conditions for a set of couplings that lead to a certain cr(H). 
Suppose that H satisfies cr(H) = c(H)\ then there must be a unitary matrix U such 


that: 


U If(N)T+ T'F > (N)\u* = F(N)T + nP i (N) 


( 12 ) 
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H = A{N)H 0 A\N), H = A(N)H 0 A\N), H 0 = T + T\ (13) 

reflecting the fact that couplings can be obtained by lateral transformations acting on 
the homogenous chain hamiltonian H 0 (if the chain is finite, the system is sometimes 
referred to as a locally periodic structure [25]). The operator A(N) is constructed by 
observing that 

A(N)H 0 A\N ) = A(N)A\N+ 1)T + T'A(N + 1)A\N), (14) 

leading to F(N) = A(N)A'(N + 1), and from here we can find the elements of A using 
the recurrence 


A n A* n+1 = F n 

or in polar form (i.e. A n = \ A n \e uin , F n = \F n \e^ n ) 

Fn | A n 11 A n +i |, f n a n o n -|_i. 

The solutions are 


(15) 

(16) 


n— 1 

^ fj, fo £ 9ft 

3=0 


(17) 


and 


Tt[|«/2|] I P I 

IA, t .h ‘S' —, i <«<JV, 

-F n 


rifc=0 \ r n-2k-L 

Furthermore, since F n can be chosen real and positive, we may work with 


(18) 


A n -y 1 


n 


[\n/2\] 


k =0 Fn-2k 
■rr[! n / 2 ll p 

Ilfc=0 “ n—2k—l 


x 4 _) ", 1 < n < N, A 1 e K\{0} (19) 


and a n = 0. Now that we have constructed A from any given F, we can also define 
its inverse restricting ourselves to the physical block, i.e. A~ l A = AA~ l = 1 ^+ 1 , with 
1jv+i the (N + 1) x (N + 1) identity matrix. With this, we can transform (12) into 


H 0 = XH 0 X\ X = A~\N)UA(N). (20) 

Our task is now simplified to the determination of all matrices X subject to (20). In 
other words, it is enough to analyze the locally periodic chain. 
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2.1.1. The matrix X The solutions of (20) are multifold. We must note that X could 
be non-unitary, which renders a solution space beyond the set of all symmetries of the 
periodic chain. We also note that the validity of (20) is restricted to the physical N + 1 
dimensional block, which is more conveniently described by N+1 dimensional projectors 
P as follows: 


XP(T + T t )P t X f = P(T + T ] )Pl (21) 

In order to find X, we compute the matrix elements of our operators in the eigenbasis 
of P(T + T')P'. For infinite systems, the states in question are Bloch waves. For finite 
systems, both wavefunctions (m\k) and eigenvalues E k of H 0 are known (in fact, they 
are related to the Dirichlet kernel). We have 


J2(k\X\m)(m\P(T + T^)P'\n)(n\X\k') = E k S k , k > (22) 

n,m 

which means that the matrix elements of X in mixed bases X k . m = (k\X\m) constitute 
an operator that diagonalizes H 0 . Therefore, we must have 


X k , m = e- ia ^(k\m), X* m]k = e ia ^(m\k). (23) 

Here we have a freedom of choice, since the phase encoded by a(k) is so far arbitrary 
§. This possibility renders X unitary, but more general choices are attainable such as 
[X, P(T + Tt)ptj = o and (XX* - 1)P(T + T+)Pt = 0, leading to 

XX t = l + P 0 , P 0 {T + T*) = (T + T t )P 0 = 0, (24) 

which is indeed more general than a unitarity condition. In fact, the projector Pq onto 
the null space of H 0 exists whenever N is even, making + 1 odd and at least one 
eigenvalue would vanish. For our purposes (general N ), the choice (23) is a reasonable 
one, as it contains a sufficient number of real parameters a(l),a(2), ...,a(N + 1). 

In addition to (21), we have the condition U = AXA~ l with U unitary. This 
implies 


(A^XA^AX^A- 1 = ljv+i 

(25) 

XA^AX* = A^A. 

(26) 


§ In the limit N —> oo, these phases are linked to the problem of the most localized Wannier functions. 
In linear arrays, this problem has been solved [21]. 
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From here we see that given a matrix A^A, we can build all possible A^A by choosing a 
specific X, which in turn amounts to a choice of the phases a{k). This expression can 
be given also in terms of matrix elements: 

|Ai| 2( 5n,n' = ^2 (n\X\k) (k\m) {m\A^ A\m')(m'\k'){k'\X^\n') 

k,k',m,m' 

= J2 e i(a(fe,)_Q:(fe)) (n\k)(k\m)(m\k')(k'\n ')| A m | 2 . (27) 

k,k',m 

The diagonal part of (27) represents indeed the orbits of a vector with components 
\A m \ 2 ,m = 1 + 1 under the action of a linear operator that depends on the 

parameters a(k). On the other hand, the off-diagonal part restricts the possible a’s as 
functions of A. This establishes that only some pairs A, A can be related according to 
restrictions such as 


y~] sin (a(k') — a(k )) (■ n\k)(k\m)(m\k , )(k , \n')\A m \ 2 = 0, (28) 

k,k',m 

| ^2 \A n \ 2 (k\n)(n\k')\ 2 = | \A m \ 2 (k\m)(m\k')\ 2 , (29) 

n m 

N+l N+l 

n ii n i 2 =n iA n i 2 , (30) 

n= 1 n= 1 

which can be easily deduced. 

To summarize this part, we underscore that the general isospectral problem for 
tight-binding hamiltonians of equal sites can be reduced to the isospectral problem of 
the (locally) periodic chain by means of a family of transformations represented by the 
matrix X. Then by varying X we obtain all the possible orbits of a given matrix A^A 
using the relation XA^AX^ = JAA = AlA[a], i.e. a functional of a. The final step 
consists in a reconstruction of F using once more F n = |A n | |7l n +i|, where additional 
phases are irrelevant due to (9). The parameter A\ is lost in the expression for F n (see 
the power (—l) n in (18)), and the generation of new couplings F n comes exclusively 
from the variation of all possible cc’s. We refer once more to figure 2. 

2.2. Algebraic equations for couplings and energies 

Here we obtain a set of algebraic equations whose solutions determine the systems that 
lead to a given (target) spectrum. Our approach in this part is a rather direct one: we 
shall find equations in the variables {F n }^Ai and parameters a(H) = {E k }k=i ■ 

We start by defining A = —E, ^(A) = det(7/ + A), leaving us with a secular 
equation <hv(A) = 0. An expansion of $jv(A) by the method of minors leads to a 
recurrence relation in the size of the system N+l: 
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3*iv(A) — A<3>7 v_i(A) — F^<3> at— 2 (A). (31) 

We want to determine the coefficients Aj v _ 1 in the following power expansion 

N+l 

= J2 A W < 32 ) 

j=0 

To this end, we replace (32) in (31) and collect terms in A- 7 to obtain 

-A-AT = ^N-l ~ ^N^N-2- (33) 

This is a double recurrence in size N and order j. In order to solve it, we need to 

specify boundary conditions on the lattice ( N,j ) (see figure 3). The index j is limited 

by N + 1, so we have A^ +1 = 1 (monic polynomial). For j = 0, we have a lower 
boundary condition provided by (33) if we set A^ v 1 _ 1 = 0, i.e. 


A 


0 _ 

N ~ 


-Fn A 


0 

N—2 


which is solved by 


(34) 


A 


0 _ 

N ~ 


0 


(TV—1)/2 

'_l)(^- 1 )/2 TT r2 


n 3 


N—2k 


k= 1 


for N even 
for N odd 


(35) 


It is important to recognize that N = Wnax hr figure 3 is not a boundary condition, 
but the mere statement that the recurrence shall be solved for a fixed (arbitrary) N. 

Now we take advantage of the fact that (33) is of order 1 in j, for the method to 
solve (33) consists in iterations of A) v _ 2 in terms A) v _ 3 , A] v _ 4 in terms A] v _ 5 , and so 
forth. This leads to a relation 

(N-j+ 1)/2 7i—l 

Mv — ( — l) n A] vr _ 2n _i Ftf_2i (36) 

n =0 1=0 

which puts A j n in terms of all the coefficients of a preceding order in j, as expected. 
Now we continue this iterative process by substituting A^ 2n _ 1 in (36) using (36) again, 
but with properly shifted indices. The process ends when we reach A^-, specified by 
(35). The formula thus obtained is 
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Figure 3. Lattice representation of the recurrence relation (33) in red and its boundary 
conditions in blue. 



with the shorthands 


0 


no=l,....,rij—i=l m =0 



if N — j even 

(37) 

if IV — j odd 


' —i— 1 ^ 

N k ,j = -|- n *’ 0 < k < j - 1, (38) 

i =0 
m— 1 

M m = N — 2 l m — 2 ni — m, 0 < m < j, (39) 

i=0 

and with n m , l m dummy indices. Summations with negative upper limits are simply zero. 
This general expression for the coefficients of the characteristic polynomial indicates that 
we are dealing with polynomials of 2iV-th order in the couplings. Also, the coefficients 
are homogeneous functions of degree N + 1 — j, and setting = constant defines a 
number of [|JV/2|] surfaces that are shape-invariant under global scaling. 

What we have achieved so far is an explicit calculation of the characteristic 
polynomial for tight-binding hamiltonians of variable couplings. Since we are interested 
in the solutions of the inverse spectral problem, we proceed now to the determination 
of <&w(A) in terms of {E k }k=x- We have 


N+l 

®n(\) = t[(X + E k ) 

k =1 


( 40 ) 
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from which the coefficients are obtained by computing the corresponding products: 

N-j 

Ajv = II Ek - M = { k > = 1 ,-,N+ lVi; hi ± kfii ± j}. (41) 

{fci} r=1 

This formula specihes uniquely the characteristic polynomial for a given spectrum; 
therefore, if we want to find the couplings F n we must equate (41) and (37). The 
resulting system reads 


(^) 


{«*,,) i 

e (-i) a; " in ' n 


hud 


m =0 



N-j 

e n < 42 > 

(fed »'= i 


This is a set of [|(N + 1)/2|] equations, giving rise to a number of N — [|(N + 1)/2|] = 
[|iV/2|] free variables and consequently a [|iV/2|]-dimensional surface representing S h- 
We shall see that these relations can be solved uniquely by providing a set of constraints 
in some of the F' s. The numerical method of choice is iterative, but we ellaborate more 
on this in section 2.5. 


2.2.1. Minimal conditions for a unique solution We have seen that specifying the 
energies leads to a surface S of models. It is possible to fix the values of [|IV/2|] couplings 
in order to obtain a single point in S . However, there is another way of defining one and 
only one model by means of the wavefunction; we show this by considering two cases. 
If N is even, the size N + 1 is odd and one of the energies vanishes. Suppose that the 
corresponding eigenvector has components with k fixed. The relation (7) implies 


F n 



iv-ir'ii 

k 


<t>n-l-k 

<t>n+l-k 


(43) 


fixing all F 1 s up to an overall scale factor F\. If At is odd, we may use our knowledge of 
any energy Ek ^ 0 to write a similar recurrence 


F n = 



(44) 


which can be solved by iteration. Once the couplings are fixed, the information of other 
wavefunctions and energies is also fixed. 


2.3. Problems with analytic solutions 

Quite often one finds that the problem of couplings for a given spectrum is accompanied 
by additional restrictions. In fact, the latter can make the inverse problem easier to solve, 
as shown in recent works [8-10]. 

The aforementioned restrictions may come in the form of Lie algebraic relations 
satisfied by the hamiltonian or by the operators contained in it. Such relations affect 
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Table 1. A class of solvable inverse problems and their algebraic structures 


Dirac oscillator 


Restrictions 


{FT, TlRt} = H 2 , 
(FT) 2 = 0 


Hamiltonian 
Group 
Spectrum 
Coupling F n 


H = cr + a + <t_<R 
S(2) 

{±V^}n S N 
\fn or 1 


Finite oscillator 

[FT, TtRt] = 2 J z , 
FT = J + 

H = J X 
SU(2) 

r TV TVl 

t 2 7 2 J _ 

\J{n - 1)(N + 1 - n) 


Position operator 

[FT, rtFt] = n, 
FT = a 


H = x 
Heisenberg 
3 ? 

\fn 


Infinite chain 


F = 1 

H = T + Tt 
(Z, +) 
[-2,2] c 3? 
1 


directly the couplings F n and impose their functional form up to gauge transformations. 
In table 1 we give some examples related to Lie algebraic structures imposed as 
restrictions. The family presented here is related to sl(2, C), which has three generators. 
This is in fact the most general structure that one may propose, since the hamiltonian 
H = FT + (FT)' is composed of only two operators FT, T^F^ belonging to a Cartan 
basis. We note however, that the possibilities provided by this well- studied object are 
sufficiently rich: compact and non-compact subsets of the special linear group give rise 
to finite or infinite spectra. 


2-4- Exact solutions for shoid chains 

A depiction of the spaces S is provided for some examples. The low dimensional cases 
illustrate in a simple manner the multiplicity of solutions for our inverse problems, which 
are presented for N = 2, 3,4, as well as the explicit constraints (42). This renders S as 
a surface or as an intersection of surfaces in the space of couplings. 



Figure 4. Curves and surfaces in the space of couplings Fi,F 2 ,..., representing 
isospectral chains, (a) N = 2, a circle; (b) N = 3, four loops at the intersections 
of a sphere with a hyperbolic cylinder; (c) N = 4, a surface generated by successive 
intersections of surfaces, each represented by a value of F 4 . 
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N = 2: a(H) = { — a JF? + F 2 2 , 0, + a/T\ 2 + invariant if F 2 + F| =constant. 

N = 3: a(H ) = {—Ei, —E 0 , +E 0 , +£’i}, with the constraints 


F? + F* + F* = E% + E%, 

F\Fl = ElE\. 

N = 4: cr(H) = {— E 2 , —E\, 0, +Fi, +F 2 }, with the constraints 

F 2 + F| + F 3 2 + F| = E\ + E 2 2 , 

F 2 F 2 + FlFl + F|F 2 = ElEl 


(45) 


(46) 


As shown before, in general, adjacent values of N (i.e. with an equal [\(N + 1)/2|]) 
give rise to sets S with the same dimensionality. It is found that N — 1 yields a circle 
and N = 2 yields a loop defined by the intersection of a sphere and a quartic hyperbolic 
cylinder. The connected components of these surfaces have the same topology. This is 
displayed in figures 4(a) and (b). When N = 4, the geometric representation is generated 
by the successive intersections of two families of surfaces in the space F \, F 2 , F :i , giving 
rise to a familiy of curves parameterized by F 4 ; as a final result a 2d surface in the space 
F \, F' 2 , F :i is obtained. This is shown in figure 4(c). 




E F 



Figure 5. Energy and coupling distributions. Upper row: (a) Cosine distribution of 
energies and (b) peaked distribution of couplings around F = 1/2 with a nonvanishing 
width. The peak represents a locally periodic system. Lower row: (c) Gaussian 
distribution for energies and (d) a distribution of couplings with a smaller variance 
and non-vanishing average. 
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2.5. Coupling distributions from random spectra: numerical examples 

When N is large, the solutions of (42) can be obtained numerically. As a demonstration, 
we generate a target spectrum from random distributions and find the family of couplings 
that produces it. This procedure gives more importance to overall trends rather than 
specific details. Our goal is to obtain the most probable couplings for two types of 
spectra: i) a dispersion relation given by a cosine law, emulating a periodic chain and ii) 
energies of a large system that are located predominantly in a gaussian band with zero 
average and adjustable width. The exploration of S is achieved by choosing random 
initial conditions in the iterative method (e.g. Newton-Raphson). 

Example i) is relevant to our general theory, given the role of X in (20). This 
offers interesting results: If the system is periodic, the spectrum must follow a cosine 
distribution law, but other possibilities for the couplings may give rise to the same 
dispersion relation, despite the lack of local periodicity. This is shown in figures 5(a) 
and (b); in (b) we find indeed a peak at F = 1/2 as the most problable value (i.e. the 
periodic case), accompanied by a nonvanishing width that reveals the existence of other 
systems with the same spectrum (note however that the eigenfunctions do not coincide). 
The average and deviation are (F) = 0.45 and (A F) = 0.20 respectively. 

The second example gives a distribution of couplings (see histograms in figures 5(c) 
and (d)) with a decay faster than exponential. This is natural, given the fact that the 
spectral span is roughly the width a of the gaussian in (c), and the maximum coupling 
of the system must be constrained also by a gaussian band, i.e. the coupling distribution 
and the energy distribution have a similar asymptotic behaviour. The average of the 
coupling distribution is (F) = 0.17 . We also find that the most probable value of the 
coupling is F* = 0.11 and it is considerably greater than zero (here 0 < F < 0.5); this is 
explained by noting that nonzero couplings may give rise to zero or close to zero energy 
eigenvalues. Thus, if the spectral distribution is dominated by the centre of a band, 
we obtain matrix elements of H that are predominantly displaced from the centre, as 
shown in figure 5(d). 

3. Bent waveguide realizations 

We test our tight-binding models with a physical system that provides the required 
flexibility. Adjustable couplings, as well as the existence of evanescent modes draw 
our attention to systems of waveguides with several bendings. Such realizations are 
indeed susceptible of manipulations that produce partial isospectrality: we should bear 
in mind that the Dirichlct boundary value problem posed by bent waveguides leads to 
nontrivial behavior at high energies or frequencies, connecting the problem with wave¬ 
like manifestations of chaos [26]. However, at very low energies we can take advantage 
of trapped states at corners [14,27-30] (in fact, only one bound state at a right bending 
angle), providing an exponentially evanescent coupling of two corners as a function of 
the distance between them. The idea then is to accomodate as many levels as necessary 
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below the propagation threshold of the guide; this is done by adjusting the distances 
between the corners contained in the array. See figures 1, 6 and 7. In this way we 
engineer the necessary couplings for a given spectrum and finally tame the beast. 

It is important to mention that, from the physical point of view, some realizations 
put restrictions on the width of waveguides [31]. Our realizations are of a different 
nature, and the existence of bound states in corners is independent of the scale. In 
fact, they have been confirmed experimentally even with microwaves [14], This can be 
easily explained by noting that the classical theorems found by Exner and Seba rely on 
the smoothness of the waveguide - in contrast to our case- as well as the independent 
variation of width and curvature. Quoting their result, “a bound state with energy 
below the Erst transversal mode exists for all sufficiently small <i”, which ensures the 
existence of a trapped mode for a certain range of the width and “the proof... establishes 
existence of the critical width do ...”. In addition, Goldstone and Jaffe (see fig. 1 in [32]) 
proved by a variational method the existence of bound states for all widths in twisted 
tubes of constant section. Right-angle bends, on the other hand, have infinite curvature 
at the corner and do not have constant section. For a single corner, one is convinced that 
increasing the separation distance between boundaries amounts to a global scaling of the 
system - see the explanation below (7) in [13]. In connection with scaling properties in 
the limit of vanishing width, some works have shown the persistence of evanescent modes 
in effectively one-dimensional wires [33], but here we provide an explicit dependence of 
the evanescence length in terms of the geometry. 

3.1. The validity of tight-binding approximations 

The equation to solve here is of the Helmholtz type (a quantum particle trapped in the 
cavity or TE modes in a good conductor): 


[V 2 + k 2 ] f(x, y) = 0, (p(x,y) = 0 if (x, y) G dVt. (47) 

The associated Hamiltonian for this equation is simply H = —V 2 . For two adjacent 
corners 1 and 2, we find the coupling F by means of an overlap between trapped modes 


F 


<f>i(x,y)H(f> 2 (x,y)dxdy t 


(48) 


giving rise to the following block representation of a two-corner system at very low 
energies: 


H] 


Low E 


E b F\ 
F E b ) 


(49) 


Here, k 2 = E b is the energy of a trapped mode in an isolated corner. Numerical 
calculations show that for an infinite guide E b = 0.93 E t [12], that is, a fraction of the 
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propagation threshold E t = n 2 /Lr if the straight segment has a width L. It has been 
shown theoretically and experimentally [14] that such trapped modes have exponentially 
vanishing tails as functions of the distance d from the corner, measured along a straight 
segment. We write 4>(x,y) ~ exp (—d/ 2A) with A ~ 1 f \fE~b- Therefore, the function F 
can be written to a very good approximation as 

F = A (d) = A(0) exp f-j 

We have verified this law numerically for a system of two corners in different 
configurations, the results are shown in figures 6 and 7. This is also in agreement 
with the classic results by Carini and co workers [30]. 

With these considerations we justify the applicability of tight-binding models for 
the description of spectra and wave functions of more intricate systems of waveguides. 
We also perform numerical computations of the 2d problem for each realization, in order 
to give a fair comparison with our predictions. 




Figure 6. Numerical results for the decay of energy splittings as a function of distance 
d in units of width L. The close agreement with the red curve is a clear indicative of 
exponentially vanishing couplings; the parameters are A = 1.3L, A(0) = 0.43-Et,. This 
has been computed for U and S configurations shown in figure 7. 


3.2. The realization of exactly solvable configurations 

Here we provide waveguide designs that produce the models mentioned in table 1, as 
well as the examples discussed in section 2.4. 

3.2.1. The Dirac oscillator This system was recently produced experimentally in [10], 
with suggestions provided in [8]. The incarnation using microwave resonators provided 
a natural way to adjust couplings by varying distances between dielectric disks. Here 
we resort to the variation of segment lengths and their orientation to build a sequence 
of dimers; the intra dimer distance has a constant value in compliance with a constant 
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Figure 7. Symmetric and antisymmetric states for S and U waveguides. The finite 
separation distance between corners produces a level splitting. 

hjijijijtjtj 

Figure 8. A bent waveguide representing a Dirac oscillator. The logarithmically 
deformed array comprises 11 dimers. Its spectrum corresponds to panel (b) of figure 
9. 



Figure 9. Spectra of waveguides representing Dirac oscillators of three different sizes: 
(a) N = 10 (b) N = 22 (c) N = 32. A better agreement with the theoretical prediction 
is achieved as the size of the array is increased. More energy levels (blue dots) can be 
fitted by a square root law (solid black curve) centered at the bound state energy (red 
dashed line). 
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Figure 10. Numerical spectra for waveguides obeying the coupling prescriptions of 
finite oscillators, (a) and (b) show results for various spectral spacings e in terms of 
C = ( E t — Ed)/l (see inset) and sizes N = 21 = 6 and 8. (c) and (d) show especial 
examples where the agreement with equispaced spectra is outstanding. 


coupling, while the inter distance of consecutive dimers is varied with a law of the type 
log y/n. For the explicit Id hamiltonian of this system, see table 1. 

The resulting waveguide and the numerically obtained spectrum are shown in figures 
9, 8. We observe an excellent agreement with our predictions: a square root law for the 
spectrum is obtained for half the energies (around the centre of the band). A parabolic 
curve is fitted to the lowest energies applying the method of least squares to the lower 
half of the spectrum||. 

3.2.2. Equispaced spectrum or finite oscillator Finite oscillators [35,36] and relativistic 
rotors [9] have been proposed in the past, with the aim of producing systems that evolve 
coherently due to their equispaced frequencies. In one of such applications it was even 
possible to produce a schematic spectrum for baryons. These constructions do not suffer 
the effects of finite size, as in the previous case. In its simplest form (see table 1), the 
spectrum corresponds to that of J x , i.e. a projection of angular momentum orthogonal 
to its quantization axis. Spectra and wavefunctions are presented in figures 10 and 11. 

|| Finite size effects can alter the upper half of the eigenvalues. For this reason we do not include them 
in the parabolic fit. 
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Figure 11. Examples of wavefunctions for a finite oscillator with 9 sites, l = 4. (a) 
m = 1, (b) m = 4, (c) m = 6, (d) m = 9. The spectrum is shown in figure 10(d). 
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Figure 12. The spectra for triads of partially isospectral systems, two sizes: (a) 
N = 3 (four sites), (b) N = 4 (five sites). Abscissa: three systems and a randomly 
generated target spectrum. Ordinate: energy eigenvalues. The arrays were obtained 
by finding the couplings and solving for the lengths of the waveguide segments. The 
boundary value problem was solved numerically and the agreement is quite reasonable. 
The uncertainties due to discretization were also considered and the corresponding 
spectrum was calculated when lengths were rounded up (dashed red lines) and rounded 
down (dashed blue lines). The isospectrality is better achieved numerically in case (c). 
The arrays are shown in fig. 13 


3.2.3. Partially isospectral systems for low N As a final application, we bnilcl a family 
of bent waveguide arrays that possess the same spectrum below propagation threshold. 
This is most easily done for low values of N, already studied in section 2.4. Three 
types of waveguides with N = 4, 5 are designed such that their spectra below threshold 
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(a) (b) 



V 


Figure 13. Wavefunctions for partially isospectral waveguides, (a), (b) and (c) 
correspond to systems 1, 2 and 3 in figure 12(b), first excited state. 



coincide. See figures 12 and 13. This example is of the utmost importance, since it 
shows that at least N + 1 levels of a waveguide (or billiard, if closed at the ends) can 
be produced by an infinite number of systems. Although the shape of the pipe cannot 
be determined uniquely using the first N + 1 levels, we have shown that the set of all 
pipes can be analytically determined for low N. 

4. Conclusion 

Inverse lattice design in one dimension shows that arbitrary symmetric spectra can be 
emulated in various types of tight binding systems. Arbitrary but finite spectra can be 
also produced by embedding it in half of the energy band obtained from a deformed 
chain; in this case a bent waveguide operating below threshold was used. The technical 
goals that we have achieved in this paper are the characterization of the space of all 
isospectral nearest-neighbour models, the explicit form of the algebraic equations that 
relate couplings and energies, as well as their analytical (short chains) and numerical 
(long chains) solutions. Moreover, we designed bent waveguides of variable segments 
that emulate a number of paradigmatic cases, such as relativistic hamiltonians (a Dirac 
oscillator) an equispaced spectrum without finite size effects (a finite oscillator) and a 
few partially isospectral 2d systems (isospectrality only below threshold). In our view, 
these results are valuable both practically and fundamentally. They can be readily used 
in order to engineer the resonances of many experiments, particularly with microwaves. 
From the mathematical point of view, we have shown the multivaluedness of solutions 
and the algorithm to obtain matrix elements from eigenvalues. 
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